Faraday waves on a viscoelastic liquid 
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We investigate Faraday waves on a viscoelastic liquid. On- 
set measurements and a nonlinear phase diagram for the se- 
lected patterns are presented. By virtue of the elasticity of 
the material a surface resonance synchronous to the external 
drive competes with the usual subharmonic Faraday instabil- 
ity. Close to the bicriticality the nonlinear wave interaction 
gives rise to a variety of novel surface states: Localised patches 
of hexagons, hexagonal superlattices, coexistence of hexagons 
and lines. Theoretical stability calculations and qualitative 
resonance arguments support the experimental observations. 

PACS: 47.50. +d 47.54.+r 47.20.Ma 

The generation of standing waves on the surface of a 
vertically vibrated Newtonian fluid (Faraday waves) is 
one of the classical hydrodynamic pattern forming insta- 
bilities [EJ. Theoretical and experimental investigations 
during the last decade have substantially improved our 
understanding of the underlying processes . Recently 
pattern formation in viscoelastic fluids came into the fo- 
cus of nonlinear science, too. The memory of viscoelastic 
fluids introduces an additional time scale, which gives rise 
to a number of interesting phenomena It is amaz- 

ing that already Faraday jjj tried to compare Newtonian 
fluids with viscoelastic ones: "The difference between oil 
and white of egg is remarkable; (...) the crispated state 
may be a useful and even important indication of the 
internal constitution of different fluids. " In Newtonian 
fluids the prevailing surface mode of the Faraday experi- 
ment is usually subharmonic (S), i.e. the waves oscillate 
with twice the period of the external excitation. Never- 
theless, for very thin layers the surface oscillation may 
become harmonic (H), i.e. synchronous to the external 
drive However, the necessary vibration conditions 

are rather extreme, rendering this parameter region dif- 
ficult to explore. 

On the other hand, viscoelastic Maxwell fluids have 
recently been predicted to exhibit the harmonic Faraday 
resonance [ fl0| , if the external drive frequency compares 
to the elastic relaxation time of the fluid. This statement 
is based on the elasticity of the material and independent 
of the filling level. Thus selecting a polymeric liquid with 
an appropriate elastic time scale will open a new way to 
a systematic investigation of the harmonic Faraday in- 
stability. The present paper deals with such an experi- 
ment. The predicted harmonic resonance competes with 
the usual subharmonic instability. Close to the bicritical- 



ity an interesting new surface dynamics occurs, unknown 
in Newtonian fluids. 

Our working fluid is a 1.5% (by weight) solution of 
PAA (polyacrylamide-co-acrylic acid of molecular weight 
5 x 10 6 , Aldrich 18,127-7) in a 60 : 40 glycerol-water 
solvent. The tabulated values for density, surface ten- 
sion and dynamic viscosity of the solvent at 20°C are 
p = 1,154/cg/m 3 , a = 0.069N/m and 77 = llmPas. 
For the polymeric solution the complex dynamic viscosity 
rf{f) = rj'(f) — irj"(f) has been measured at room tem- 
perature (~ 22°C) with a rotational viscometer (Rhco- 
metrics Fluids). Data acquisition for frequencies above 
/ = 16Hz was not possible. However, within the range 
1Hz < f < lQHz the data could well be fitted by 
the following power laws: rf = 0.082 x f~ 594 Pas and 
77" = 0.122 x f~ 549 Pas, where / is measured in Hz. 
These relations have been used to extrapolate the vis- 
cosity data into the frequency range of our experiment 
(30Hz < f < 100Hz). 

The experimental setup consist of a black container 
built out of anodised aluminium with an inner diameter 
of 290 mm, filled to a height of 3 mm with the work- 
ing fluid. By means of an electromagnetic shaker (V617 
Gearing & Watson) the container is vibrated vertically 
with an acceleration a cos fit, where f2 = 2irf. The con- 
tainer is covered and sealed by a glass plate to prevent 
evaporation and pollution. No degradation or drift of the 
experimental results could be detected within a period of 
four weeks. Due to the heat production of the vibrator 
we were not able to work at room temperature, where 
the rheometric viscosity data have been taken. A heat 
wire filament regulates the temperature at 30°C ±0.1, 
measured by a PTC-resistor imbeded in the bottom of 
the vessel. Waveform generation as well as data acqui- 
sition are performed by a D/A-board (National Instru- 
ments, AT-MIO 16 E2) in a Pentium PC. The vibration 
signal measured by a piezoelectric device (Bruel & Kjear 
4393) was checked to be sinusoidal with contributions 
from higher harmonics not more than 2%. A ring of 150 
LED's with twice the container diameter is mounted on 
the vertical axis of the setup and illuminates the surface. 
The light reflected into the centre of the ring is recorded 
by a full frame CCD camera (Hitachi KPF-1). Its elec- 
tronic shutter allows to synchronise the exposure with the 
external drive signal an thus to discriminate the H and S 
surface response. The phase diagram of the selected pat- 
terns is obtained by quasi-statically ramping the drive 
amplitude from just below the onset e = ^ — 1 = —5% 
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FIG. 1. The acceleration amplitude a c and critical wave 
numbers k,H and k$ for the onset of the harmonic (circles) 
and subharmonic (squares) Faraday instability. Experimen- 
tal data are plotted by symbols, lines denote the theoreti- 
cal computation on the basis of the rheometric viscosity data 

»?*(/)• 

up to e = 10% while keeping the frequency constant. 
During each scan the amplitude is increased by 1%-steps 
and held constant for 300 seconds in between. Then, 
the surface of the fluid is photographed and the time 
dependence of the surface oscillation (whether S or H) 
is determined. Similar downwards-amplitude scans are 
performed to detect a possible hysteresis. 

Fig. [l] shows the threshold amplitude a c for the critical 
wave numbers of the Faraday instability. In the frequency 
domain above / ~ 60Hz one observes the usual subhar- 
monic Faraday resonance with a surface pattern of lines. 
Since this behavior is similar to high viscosity Newtonian 
liquids QJl^], we denote this region as the Newtonian 
regime. The present paper, however, is focused on the 
frequency range below AQHz, where the surface responds 
harmonically. The predominance of the harmonic Fara- 
day resonance has recently been predicted for viscoelas- 
tic Maxwell fluids |l(]] . On using the rheometric viscosity 
data for our working fluid we have computed the stabil- 
ity threshold (lines in Fig. by the method of Kumar 
and Tuckerman Favourable agreement is achieved 

except for frequencies close to the bicritical intersection 
point. There are two possible sources for discrepancies: 
(i) There is a 7 — 9°C difference between the tempera- 
ture of our experimental setup and that of the rheometer 
during the viscosity measurement, (ii) Errors due to the 
extrapolation of the viscometric data into the frequency 
range of our experiment. 

Fig. ^ depicts the phase diagram of the selected pat- 
terns. The surface structure in region Ila is a uniform 
pattern of hexagons covering the whole surface, while in 
lib we observe very stable localised patches of hexagons 
(see Fig. ||a). A slight onset hysteresis up to 6% in e is 
visible. The pattern selection processes for 
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FIG. 2. Phase diagram of the observed nonlinear patterns. 
The symbols mark experimental data points. The abscissa, 
e = 0, indicates the linear threshold shown in Fig. |l| The 
lowest solid line denotes the saddle point of the hysteresis. 
Region I: flat surface; Ila: harmonic hexagons covering the 
whole surface; Ilia: harmonic-subharmonic hexagonal super- 
lattice (see text) extending over the whole surface; lib, 1 1 lb: 
as before but localised patches surrounded by the flat sur- 
face; IVa, IVb, IVc: chaotic dynamics of subharmonic lines 
competing with extended (a) or localised (b) hexagonal su- 
perlattices or the flat surface (c); V: stationary subharmonic 
lines extending over the whole surface (Newtonian regime). 

the backwards bifurcating hexagons in the harmonic 
region 77 and the forwards bifurcating lines in the sub- 
harmonic Newtonian regime V rely on the distinct re- 
spective time dependencies. This can be seen as follows: 
On taking the surface elevation r\ as a representative or- 
der parameter one has 
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Here r = (x, y) is the horizontal coordinate, the lat- 
eral wave vectors k# m and ks m with |kjj m | = kn and 
ks m | = ks compose the spatial pattern and the r/ n are 
the temporal Fourier coefficients determined by the lin- 
ear stability problem (components of the Floquet eigen- 
vector). A similar ansatz holds for the velocity field v. 
On feeding n and v into an arbitrary quadratic nonlin- 
earity of the hydrodynamic equations results in a fre- 
quency spectrum of integral multiples of f2, no matter 
whether S or H is considered. Thus quadratic nonlin- 
earities are able to resonate with harmonic linear eigen- 
modes, but not with subharmonic ones. Clearly, spatial 
resonance must be granted as well. Thus any triple of 
harmonic modes {ksi, kf/2, ^113} with |kff m | = ku and 
k#i +kf/2 + k#3 = is in resonance. This generic 3- wave 
vector coupling is well known e.g. from Non-Boussinesq- 
Rayleigh-Benard convection and enforces a saddle node 
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FIG. 3. Localised stationary surface patterns of harmonic 
hexagons (a) and the harmonic-subharmonic hexagonal su- 
perlattice (b) as observed in region lib and 1 1 lb at / = 37Hz 
and e = —0.05 and —0.03. The corners of the picture show 
the outer edge of the container. In (c) and (d) the respective 
spatial Fourier spectra are shown. Patterns in regions Ila 
and Ilia are similar but they cover the whole surface. 



bifurcation towards hexagonal patterns. The associ- 
ated amplitude equations are of the form (use cyclic per- 
mutation to get the other ones; n and 7 arc nonlinear 
coefficients) 

d t H x = eH x + kH*H* 3 - {Iff! | 2 + 7 [\H 2 \ 2 + \H 3 \ 2 ] } H x . 

(2) 

Within the Newtonian region V such a resonant 3-wave 
vector coupling is prohibited due to the missing tempo- 
ral resonance. Quadratic nonlinearities in the evolution 
equations for the S m are prohibited and the pattern se- 
lection mechanism is controlled by the cubic nonlinearity 
|lj,^2|. It is instructive to point out that the subhar- 
monic temporal symmetry r)(t + 27r/51) = —T)(t) for Fara- 
day waves is analogous to the spatial up-down symmetry 
for Rayleigh-Benard convection (known as the Boussi- 
nesq symmetry): They both prevent quadratic nonlin- 
earities to appear in the associated amplitude equations. 

If the drive amplitude is raised from Ila, b to Ilia, b, 
a sharp secondary transition can be detected: Subhar- 
monic frequency contributions suddenly appear in the 
temporal Fourier spectrum of the surface elevation. Si- 
multaneously the associated spatial Fourier spectrum 
(Fig. |^c,d) exhibits a hexagonal superstructure with the 
same orientation but a wavelength twice as long as the 
primary one. The theoretical stability calculation indi- 
cates that the primary subharmonic onset (dashed line 
in Fig. |l|) is situated far above the 77 — III transi- 
tion. Therefore, the secondary superstructure has to re- 
sult from a nonlinear excitation process. The underly- 
ing mechanism is again a 3-wave vector interaction: In- 
voking basic resonance arguments the subharmonic mode 
ksi = ilcai with amplitude Si obeys the following am- 
plitude equation (use again cyclic permutation to get the 
equations for 5*2 and S3) 



d t Si = —XSi + xHiS{* + higher order terms. 



(3) 



Here A > reflects the linear damping and x is a non- 
linear coupling coefficient. According to Eq. (||) the sec- 
ondary crossover from II to III occurs when the primary 
pattern amplitude \H m \ = H exceeds the threshold A/|x|- 
The dashed lines in Fig. || separating the subregions of 
homogeneous (a ) and localised (b) hexagons are to be 
considered as approximate boundaries rather than sharp 
transition lines. On increasing /, the localised patches 
decrease in size until they dissappear at / > 39Hz. How- 
ever, on keeping / fixed their size is rather robust and 
hardly depends on the forcing amplitude a. Note that 
these new localised patches of hexagons cannot be ex- 
plained by a set of Ginsburg-Landau equations, which 
supplements Eqs. (Q) by diffusive spatial derivatives: 
Within this familiar model isolated islands of hexagons 
do not stably exist over a finite control parameter inter- 
val. 
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FIG. 4. Snapshot of a time dependent surface state in 
region IVb at / = 37Hz and e — 2. A stationary localised 
patch of hexagons is surrounded and competes with subhar- 
monic lines moving in an erratic manner across the surface. 

On entering region IV a or IVb the patterns become 
chaotically time dependent. Patches of subharmonically 
oscillating lines originating in an erratic manner from the 
cell boundary (a) or respectively from the flat surface 
(b) penetrate into the stationary hexagonal superlattice. 
Then they disappear again and the original structure is 
recovered. This dynamics repeats on time scales of the 
order of seconds to minutes, leading to a temporary co- 
existence of the stationary hexagonal superlattice with 
subharmonic lines (see Fig. ^). Higher drive amplitudes 
lead to a fully chaotic surface pattern. 

For drive frequencies in the intermediate region Vc a 
similar chaotic behavior of subharmonic lines can be ob- 
served, however directly as the primary instability. Since 
hexagons do not occur beyond / ~ 39Hz, the subhar- 
monic lines compete with the flat surface state. Never- 
theless, the bifurcation still exhibits a small hysteresis, 
which continuously goes down as / approaches 60Hz. 

In summary, we have presented experimental results on 
Faraday waves in a viscoelastic medium. Earlier theoret- 
ical predictions of the harmonic Faraday resonance could 
be confirmed. Our measurements are the first system- 
atic study of harmonic Faraday waves, which could not 
be performed yet for Newtonian fluids. The empiric onset 
data are found in good quantitative agreement with the 
theory. Furthermore the appearance of uniform hexag- 
onal patterns and their (secondary) transition towards 
a harmonic-subharmonic superlattice is understood in 
terms of elementary 3— wave vector resonances. Besides 
spatially homogeneous wave patterns we find very robust 
and sharply localised patches, which are surrounded by 
the flat surface or - at higher drive amplitudes - in com- 



plicated dynamical competition with subharmonic lines. 
Their explanation requires a more elaborate theory. 
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